Pfaffian State Generation by Strong 3-Body Dissipation 
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We propose a scheme for preparing and stabilizing the Pfaffian state with high fidelity in rapidly rotating 2D 
traps containing a small number of bosons. The goal is achieved by strongly increasing 3-body loss processes, 
which suppress superpositions of three particles while permitting pairing. This filtering mechanism gives rise to 
reasonably small losses if the system is initialized with the right angular momentum. We discuss some methods 
for tuning 3-body interactions independently of 2-body collisions. 
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Trapped cold atoms represent a new frontier for the study of 
many body phases, due to higher versatility if compared to or- 
dinary condensed matter systems [1]. Confining an atomic gas 
to live in a 2D geometry represents a benchmark for testing 
quantum Hall effects (QHE), where fast rotation for bosons 
plays formally the same role as the magnetic field for electrons 
confined in eteroj unctions [2]. Strictly speaking, the analogy 
with the QHE is attained when the centrifugal force equals 
the radial trapping force, recovering the big degeneracy of the 
lowest Landau level (LLL). Nonetheless, some phases can be 
stabilized keeping the rotation below this limit y|]. 

One of the most intriguing states appearing in the QHE is 
the Pfaffian, first introduced by Moore and Read [4] in the 
context of paired Hall states. Successively, it was proposed 
as a candidate for describing the filling factor v — 5/2 in 
fermionc QHE, although the question is still under debate and 
the Pfaffian state has not been observed so far. However, this 
special wavefunction captures the attention of the scientific 
community, since it has the peculiar property of non-Abelian 
braiding statistics of its anyonic excitations. This feature dis- 
criminates the Pfaffian from other QHE states (like the Laugh- 
lin) and makes it attractive in the context of topological quan- 
tum computation (TQC) lUS]. 

The aim of this Letter is to propose and investigate a novel 
method for preparing and stabilizing the Pfaffian state in rotat- 
ing harmonic traps loaded with cold bosonic atoms. Our pro- 
posal is based on exploiting dissipation instead of suffering 
from it 121 l8l |9[], by implementing a filtering procedure that 
basically projects a state with the right angular momentum 
onto the desired Pfaffian. In experiments, this can be achieved 
by a readily feasible mechanism to increment the relative im- 
portance of 3-body losses with respect to 2-body elastic scat- 
tering: namely, squeezing the trap lioll . Furthermore, perfect 
projection is attained by turning off 2-body collisions (by Fes- 
hbach resonances 111 1I0 . The method of preparing the initial 
state with a given angular momentum is not crucial for the 
filtering to work, but we expect it to be possible in the near fu- 
ture by putting few atoms in single wells of an optical lattice 
into rotation [12], or by realizing gauge potentials [13]. 

The paper is organized as follows. First, we briefly re- 
vise the QHE regime requirements for cold bosons and the 
physics of the LLL, which is dominated by the nature of in- 
teractions. Hereby we demonstrate the different behavior of 



2- and 3-body contact potential under density rescaling in 2D. 
The ground state (GS) diagram in the case of conservative 3- 
body scattering is presented, with particular attention to the 
Pfaffian-like sector. This motivates the need for an effective 

3- body repulsion. The dissipative projection mechanism is 
then proven without 2-body collisions, and investigated nu- 
merically in their presence. Finally we suggest some experi- 
mental signatures and a method to prepare the initial state. 

At low energy, a realization of an effective 2D system is 
obtained by setting the longitudinal trap frequency much big- 
ger than the transverse one, co± ^> CO, in order to freeze lon- 
gitudinal motion in the GS. The filling factor is defined as 
V = N/l m ax with / max the maximum angular momentum occu- 
pied by single particles. In the frame rotating at angular speed 
Q.Z, the single body Hamiltonian in the trap can be written as 



trap 



2m 



. + -(co z -Q. z )(x I +y z ), 



with A = mQz x r. In the limit of centrifugal deconfinement 
i2 — > CO, only the Coriolis force is remaining and the system 
is formally equivalent to bosons of charge q in uniform mag- 
netic field B — (2m£l/q)z. The one body eigenfunctions in the 
LLL take a simple form when written in terms of the complex 
coordinate z = (x + iy) / % , with S, = y/W/mco: 
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and have energies E„ = Tin(co — Q.) = l„ 8co where l„ are the 
angular momenta projections. Other Landau levels (LL) are 
separated by a gap ~ 2co, and the LLL restriction is claimed 
to be even more valid for higher-L QHE states, due to their 
stronger correlations [ 14]. All the numerical calculations pre- 
sented below have been checked by including the first LL and 
verifying that occupation is small there. 

The physics in the LLL is dominated by the nature and the 
strength of interactions, which drive the system into different 
filling factors |0]. In the context of cold bosonic gases in 2D, 
2-particles interactions can be modeled by contact potentials 



(1) 



with g2 



2D 



STthco^ fl/^x, being a the s-wave scattering 



length in 3D and £,± = y / h/mco± the longitudinal trap size 



2 



(]]]. Analogously, pointlike 3-body interactions can be written 
in the following way lfl5l[l6ll 

^i=gf L 8®{jtt-jtj)8&>(jtj-jt k ), (2) 
i<j<k 

neglecting for the moment their physical origin. Let us call K2 
and K3 the kernel of 2-body Eq.(Q]) and 3-body Eq.© term, 
respectively. Inside K2, the Laughlin state with V = 1/2 



^Lau-nfe-^') 2 : 



(3) 



has the lowest total angular momentum LLau —N(N—1), and 
thus energy £Lau = ^LauSo. [As usual, in © and in the subse- 
quent QHE wavefunctions we omit the ubiquitous exponential 
and normalization factors]. Of course, K2cK3 and (O is also 
annihilated by the 3-body interaction, despite not being the 
GS. The lowest L state in K3 is indeed the Pfaffian 



«Ppf = Pf 



Zi-Zi 



(4) 



withLpf = N(N-2)/2 and v = 1. The pref actor Pffjpjr) 
makes possible the superposition of pairs, and it is formally 
equivalent to a projected p + ip BCS wavefunction jl6lll7ll . 

The interplay between 2- and 3-body terms and the imper- 
fect matching 8 CO > of trapping and rotation frequencies 
enriches the phase diagram, by stabilizing other states aside 
(OH 1 3]. The sequence of GSs has been computed numer- 
ically by exact diagonalization in the case of elastic 3-body 
repulsion; the result is presented in FigQ]for N = 6, in a trun- 
cated LLL single particle basis that contains all the angular 
momenta up to / max = 2(N — 1 ) . This is enough to describe ex- 
acly ^Lau* hence suitable also for the description of states with 
lower angular momenta. In the LLL, central contact potentials 
(Q]|2]) can be expressed in terms of a single pseudopotential 
acting only on zero relative angular momenta 11811 . Thus their 
sum can be recast in a very compact form lisill 



(5) 
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where d„i is an annihilation operator of n particles with total 
angular momentum projection /. The GS energy Eq(L) is a 
monotone nonincreasing function of L. The global GS simply 
belongs to the sector of L that minimizes the quantity Eq(L) + 
L8co. Notice that the separation lines in Fig. Q]are nearly 

straight. In absence of the 3-body term, C3 = 0, the GS of J#2 

(2) 

with L = Lpf = 12 - that we indicate as \ Yq (12)) for brevity 
- is unique and stable for a narrow interval of 8(0. The state 

\Yq (12)) and the Pfaffian state (0]i share some similarities, 
since they have the same angular momentum and their relative 

fidelity isF= |( v P Pf |v/d 2) (12))| 2 ~0.803. The most important 
feature of 3-body interaction is the enlargement of both the 
stability interval and the fidelity, up to F ~ 0.991, for C3 = 1. 

The previous analysis urges us to design a mechanism for 
independent tuning of 2- and 3-body terms, a task we tackle 



,3 0.3 




Figure 1 : GS phase diagram for N = 6 particles and ci = 1 /6. Differ- 
ent colors correspond to different angular momenta of the global GS. 
We identify the Laughlin state (v = 1 /2) with L = N(N - 1) = 30, 
the Pfaffian (v = 1) with L = N(N - 2)/2 = 12, the single vortex 
with L = N and finally L = 0. 



with the following dimensional argument. In accordance 
with the fact that 5' 2 ' has dimensions of a inverse squared 
length [£~ 2 ], the coupling g 20 in (fl]i has the dimensions of 
an energy times a squared length. So, under rescaling E — > 
E/h(0, i — > £/% the adimensional coupling 02 — VSna/^j_ 
does not depend anymore on the planar trap frequency. In 
other words, for 2-body collisions 2D are special since den- 
sity does not discriminate between weakly and strongly in- 
teracting regime. Along the same line, the adimensional 3- 
body coupling scales as the inverse of the effective trap area, 
i.e. c 3 = gl°i;- 4 /hco = {gfm/h 2 )t,' 2 . This suggests that the 
relative importance of 3-body collisions can be boosted by a 
squeeze of the 2D trap, which increases the density. 

Unfortunately, elastic collisions involving 3-body processes 
are rather rare in nature. The most prominent 3-body collision 
process known in physics of bosonic condensates is due to 
recombination 12011 . The formation of a biatomic molecule 
is assisted by a third particle that assures energy-momentum 
conservation. The recombination rate displays a rich behav- 
ior as a function of the 2-body scattering length a, typically 
tuned via Feshbach resonances. For a > it shows a universal 
behaviour ~ a 4 [21], whereas for a < some genuine 3-body 



resonances are appearing due to Efimov trimer states 12011 . In- 
terestingly, 3-body processes are still present even in absence 
of 2-body scattering (a = 0). This nonlinear features allows 
to tune C3/C2, in conjuction with the squeezing mechanism 
exposed above. Tipically, recombination is considered an un- 
wanted effect in experiments with condensates since it yields 
to severe 3-body losses. Nonetheless, strong dissipation has 
been exploited successfully to induce strong 2-body correla- 
tions in the ID Tonks-Girardeau gas, as observed in recent 
experiments 10, Ht]. Moreover, 3-body dissipation has been 
proposed for obtaining a dimer superfluid phase in ID attrac- 
tive bose-Hubbard models [9]. On the same line, we propose 
here to use strong 3-body recombination rate to filter out the 
Pfaffian state (|4j, with high fidelity and paying the small price 



3 



of moderate losses. 

The Markovian dynamics of the system is described by a 
Lindblad master equation for the density matrix p 



Y^Y^dvpd],, (6) 



where 73/ is the rate of decay in the channel of 3-body total an- 
gular momentum /. The effective non-hermitian Hamiltonian 
is given by 

J^eft — ^trap + ^2 — 2 52 73/^3/^3/ = ^trap + ^2 ~ ^ ^3 • 

(7) 

In the quantum jump approach 11221 12311 . loss events are as- 
sumed to be rare and are ideally detected by performing very 
frequent measurements. This allows to describe the non- 
unitary evolution of the system in terms of wave functions, 
through the equation iti \ y/(t)) — 1 y/(t)) . The reduction of 
|| y{t) 1| 2 gives the probability of having suffered from a jump, 
i.e. a 3-body loss, in [0,t]. Actually, we are interested only in 
lossless samples and we discard the cases where losses occur, 
since generally they lead to projection onto excited states. Ex- 
perimentally, the discarding procedure is made possible by us- 
ing post-selection of the samples, after measuring the number 
of particles. In the limit case c? = 0, the dynamics is governed 
only by the dissipative term and |y/(f)) = exp(-tJ^) |y/(0)). 
Due to the positiveness of M^, this evolution projects onto K3 
for long enough times, whose scale is set by 1 /C3. 

Let us assume now to have prepared several copies of a few 
body system with predominant 2-body collisions in the GS at 
angular momentum Lpf. We will discuss later the experimen- 
tal feasibility of this assumption. At t = 0, strong 3-body dis- 
sipation is suddenly switched on by squeezing the trap. The 
ideal situation is accomplished by turning off completely c 2 
while having C3 not too small. In this case the filtering pro- 
duces the Pfaffian state since it is unique in K3 for Lpf . This 
works provided the starting \y/^\Lpf)) has a sizable overlap 
with the Pfaffian, which indeed occurs for a moderate N. 

Possibly, a more feasible experimental procedure would 
avoid to switch off c 2 at t = 0. The presence of 2-body colli- 
sions contrasts the formation of pairing contained in the Pfaf- 
fian. However, we expect that for C3/C2 3> 1 the perfect pro- 
jection is almost recovered. In order to check this statement, 
we have numerically simulated the evolution for N = 6 by us- 
ing a fourth-order Runge-Kutta algorithm. In the left panel of 
Figf2]the population of trajectories unaffected by jumps and 
the fidelity with the exact Pfaffian state (0|i are plotted as a 
function of time for different values of C3. It emerges the 
following scenario: strong enough dissipative rate yields to 
intense losses up to a time after which the population is sub- 
stantially stationary. The stabilization of losses is a signal that 
the filtering procedure has converged to a state very close to 
the Pfaffian one, as witnessed by the saturation of fidelity to 
a value close to 1. This behaviour may be interpreted as a 
sort of Quantum Zeno effect, where strong dissipation freezes 
the system in K3, suppressing losses [7]. The typical time we 
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Figure 2: Left panel: fidelity with the Pfaffian and surviving pop- 
ulation as a function of time, after switching on dissipation C3 on 

the initial state | v^ 2) (12)). The fidelity reaches 0.9994 for c 3 = 10. 
Right panel: threshold time for having lost 25% of population after 
switching the dissipation C3. 



have to wait for reaching the threshold of 25% losses is shown 
in the right panel of FigfJ] As expected, losses are not com- 
pletely suppressed for long times, since the Pfaffian is not an 
eigenstate of J%2. 

The considerable improvement in the reproduction of the 
Pfaffian state, with infidelity sinking from 20% to almost 
0.05%, though relevant by itself, does not exhaust all the im- 
portance of the proposed scheme. Indeed the flurry of interest 
about this state is mainly related to its zero-energy excitations, 
that have non-Abelian braiding properties and may constitute 
the basic ingredients for TQC [6]. We now show that our fil- 
tering procedure allows for the production and manipulation 
of the typical "half-flux" quasiholes 

l Ppf+2qholes = Pf [ [(Zi - Zj) , 

V Z i Z j / i<j 

(8) 

otherwise inaccessible with only rotation and conservative 2- 
body interactions. The quasihole identification and motion is 
indeed possible only in presence of an appropriate gap pro- 
tected subspace, i.e. states in K3 with L > Lpf. Such states 
are obviously steady states of the 3-body dissipation, and the 
quasi-constant gap guarantees that the filtering have almost 
the same speed and neatness for all of them. Then, starting 
from having filtered the Pfaffian out of an initial state, it would 
be in principle possible to engineer some excitation scheme 
being sure to lie inside K3. However, due to the 80% similar- 
ity between the GS \\]/Q 2 \L Pf )) of J^ 2 and the Pfaffian |v%), 
one might be tempted to speculate an approximate scheme for 
quasiholes with a similar precision. This is not the case, since 
the degeneracy of the quasihole subpace is completely spoiled 
out, to the point we cannot speak about a manifold protected 
by a gap. For testing our assertion, we choose the state ([8]) 
with quasiholes located at the center, i.e. w\ = W2 = 0, liv- 
ing in the sector L = Lpf + N. By expanding ^ in terms of 
the eigenvectors of J%2, we discover that this peculiar state 
has a sizable overlap with several excited states of J%2, whose 
energies are spread over an interval of the order of 
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Figure 3: Map of the gap between GS and first excited state as a 
function of the trap parameters (8(0, e). The white arrow indicate 
the chosen path for conneting adiabatically the GS with L = 6 and the 
GS with L = 12. The the matrix elements |(i)^)|ff[Vfe)| ( as needed 
by the adiabatic criterion) give a very similar picture. 



For completeness, we want to discuss about the experimen- 
tal detection of the Pfaffian state, that naturally follows the 
preparation. Of course, a great evidence is provided by the 
suppression of losses, that signals the absence of local 3-body 
superpositions. The probe is conclusive after a measurement 
of L and N, since the unicity of the Pfaffian in K3. The Pfaf- 
fian is also expected to have an increased pairing with respect 

(2) 

to iVo (12)). As a matter of fact, the expectation value of 
the 2-body contact potential is respectively 3.13 C2 and 
3.31 C2, with a difference of about 5%, maybe not striking 
enough to be resolved in a clean way in experiments. 

In the previous study, we assume that before applying the 
filtering via dissipation, it is possible to initialize the system 

(2) 

in |i//q (Lpf)). To obtain a given filling factor v it is nec- 
essary to apply a rotation to the ultracold gas. The widely 
used experimental technique consists on stirring the conden- 
sate 124 12511 . In the rotating frame, it is equivalent to induce a 



small quadrupole deformation to the trap H — > H + H £ , with 



e(x 2 -y 2 ), 



(9) 



where £ is meant to be small for avoiding heating and coupling 
with higher LLs. The term (O couples single particle states 
with angular momenta differing by 2. A possible route now 

is to find an adiabatic path in the plane (8(0, e) connecting an 

(2) 

already realized state, such as the sing le vortex GS \ % '(N)), 

12) 

to the GS in the desired sector (Lpf)). Once again we 
resort to the case of N = 6 and diagonalize the problem in the 
full LLL Hilbert space in order to draw the map of the first 
excitation gap 01211 . as displayed in Figj3] A good path would 
try to avoid regions with a small gap, since they are roughly 
associated to a slow down of the adiabatic evolution. We have 
studied numerically the time evolution along the path marked 
by the white arrows in Figj3] with an adaptive method that 
adjusts the parameter speeds according to the adiabatic con- 



dition |(y/b|H|y/ m )| < \E m -Eo\ z , Vm 7^ 0. At the end of the 
path, it is possible to reach a fidelity of 99.7% with the GS at 
L = 12. The price to pay is a evolution time t ss 10 3 © -1 which 
is quite long, but still comparable with the typical time of such 
experiments in traps. The reason is that for reaching high an- 
gular momenta, it seems unavoidable to cross regions with a 
small gap. The way how Figj3]scales with N is not clear, and it 
may be the case that the mean-field and QHE regimes are sep- 
arated by a quantum phase transition. This fact could explain 
why so far experiments failed in reaching the QHE regime 
lll2l l24l 12511 . However, for a small number of particles like 



in our case, the gap is still sizable and we expect experiments 
with the single holes of optical lattices 111 211 to succeed in the 
next future. Any other scheme of GS preparation would be 
equally good to provide the starting point for our procedure. 
It is remarkable that the filtering strategy is working best if 

(2) 

C3 is turned on after the preparation of | y/g (Lpf )) ■ Our simu- 
lations indicate that switching on C3 slowly or along the path 
for preparing \\]/^\Lpf)) produce more losses. Finally, we 
emphasize that starting the adiabatic preparation from states 
with higher L, like ^Lau, needs evolution times which are too 
long, due to crossing of regions with very small gaps. 

In conclusion, we propose a feasible filtering scheme for 
the realization and stabilization of Pfaffian state and excita- 
tions. Once the present considerable efforts about preparing a 
system in the desired angular momentum would achieve their 
goal, any other requirement is far within the present technolo- 
gies. We showed that a tuning mechanism for increasing the 
relative importance of 3-body losses relies indeed on squeez- 
ing the harmonic trap. Experiments in this direction have been 
done recently ifioll . Reducing the magnitude of 2-body colli- 
sions, by Feshbach resonance techniques, further enhances the 
filtering procedure. The almost perfect cancellation of ci have 
been obtained in Ref.[ll]. Detailed calculations and discus- 
sions are included in a forthcoming work 11911 . 
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